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ABSTRACT 

This report is concerned with the development of a new approach 
to the problem of stability for multidimensional, causal, recursive, 

*all pole*, digital filters. The distinguishing feature of this approach 
is that general stability criteria can be derived directly in terms of 
the coefficients of the transfer function of the filter. Thus by use of 
this method it is sometimes possible to determine which coefficients 
of the transfer function are critical to the stability of the filter, 
information which is, of course, important in filter design. Also the 
emphasis of this approach is on the development of a conceptual method 
for considering the problem in complete generality. 

Reproduction of all or part of this report is authorized. 
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1. Introduction 

Consider a filter whose transfer function has the form 

H(z) = 1/(1 - Q(z‘b) (1.1) 

where z = z^) is a vector complex variable, and QCz) is a 

polynomial in N variables with real coefficients and zero constant term. 
Such a transfer function describes an ’all-pole’ ’causal’, recursive, 
multidimensional digital filter, and it is known that questions of stability 
generally reduce to a question of stability for filters of this type. It 
is also known that such a filter is stable iff 

P(z) = 1 - Q(z) (1.2) 

has no zeros S = •••> S^) such that |S^| ^1 for all i . For 

this reason, it is convenient terminology to state that a polynomial P(z) 
of several variables is stable if P(S) = 0 implies that |b^| > 1 for 
some i , The question of stability then becomes, primarily, the problem 
of determining if such a polynomial is stable. 

An N-tuple of the form jf = (j^, ..., where each is a non- 

negative integer is called a vector index . If for such an index X > 
define ‘z’^ by 




then every polynomial P(z) in N-variables can be written in the form 



1 



(1.4) 



P(z) = ^ ^ 

-+ 

J 

where for each ^ , is a real numberj equal to zero except for a 

finite number of indices (see Davis and Souchon [1975]). For example if 

2 

P(z^,Z2) = 1.00 + O.Sz^z^ - .llz^z^ 



then a^Q = 1.00, a^^^^ = .05, a 2 -j^ = -.11 and all other coefficients equal 

zero. 

Suppose that a fixed, ordered set (j*. , ..., "j" ) of non-zero, 

± Z K 

N-dimensional vector indices is given, and consider the set of polynomials 
of the fonn 



P(^ = 1 




(1.5) 



where, for each i , a^ is a real number. Each such polynomial (1.5) 

■^i 

is naturally associated to the point (au^ , . . . , a^ ) of K-dimensional 

^1 

euclidean space. Conversely, each point (A^^, •••, 
euclidean space is associated to a polynomial of type (1.5) by the 
correspondence 



A^) of K-dimensional 



A « a ^ j x~l,..«,K (1*6) 

relative to the ordered set of indices . • • , • Thus relative to 

i K 

this ordered set of indices, K-dimensional space becomes the coefficient 
space for polynomials of type (1.5). A point (A^ 

coefficient space will be said to be a stable point if its associated 
p o lynomial , namely 






of this 



2 



( 1 . 7 ) 



p(2) = 1 - A^. ~z 

i=i 

is stable. The set of all such stable points in the coefficient space will 
be called the region of stability . The problem that will be considered 
here is what can be found concerning the region of stability for a given 
set of vector indices, or equivalently, for a given type of polynomial. 

In certain cases, the region of stability can be specified completely. 
Consider the ordered set of indices = (10), 'j *'2 = (01) “ (H)* The 

associated polynomials are of the form 

P(z^,Z2) = 1 - - k^z^ - k^z^z^ , (1.8) 

and the associated coefficient space is 3 dimensional euclidean space. 

Huang [1972] has shown that (1.8) is stable iff 

IA3I < 1 

IA1+A2I < 1- A3 

|a^ - A2I < 1 + A3 . 

The region of stability is illustrated in Figure 1. Also, it can be 
shown, see Jury [1974] that a similar type of specification can be given 
for any region of stability associated to polynomials in one variable, or 
equivalently, to any ordered set of 1-dimensional vector indices. 

Although, regions of stability can be very difficult to determine for 
more complicated types of multivariable polynomials, in the following it 
will be shown that several properties of regions of stability can be 
derived in general. 
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As 




FIGURE 1 



4 



2 . 



The Basic Theorems 



Throughout the following it will be assumed that "J. , . . . , ^ is a 

X K 

given, fixed, ordered set of N dimensional vector Indices. Thus the 
correspondence between polynomials of the form (1.7) and points 
(Aj^, . . . , Aj,) of K-dimensional space is assumed fixed. 

The order of a vector index "f = (j^, ..., j^) is defined by 

order j = j ^^ + . . . + . 

A coefficient A^ of a polynomial of type (1.7) is called a leading 
coefficient if its associated vector index "t. has maximal order. A 

-’i 

polynomial of several variables may have several leading coefficients. 

Theorem 2.1 If (A^, . . . , A^,) is a stable point of the coefficient space 
and if A is the sum of the leading coefficients of the polynomial 
corresponding to this point, then |a| < 1. 

Proof. Let PC?) be the stable polynomial associated to (A^, ..., 

It follows that the one variable polynomial p(z) defined by 

p(z) = P(z, ..., z) (2.1) 

is also stable. Moreover the leading coefficient of p(z) is seen to be 
A . The polynomial p(z) factors, 

p(z) = A(z-t^) . . . (z-t^) (2.2) 

where t^, ..., t^^j are its complex roots. Since p(z) is stable it 
follows that 

|t^| >1 i = 1, ..., M . (2.3) 
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But then. 



1 = |p(o)| = |a| 



M 

I n 



i=l 




(2.4) 



where 



M 

I n t I > 1 . 
i=l ^ 

Hence, 



1a| < 1 . 

Theorem 2.2 (Necessity Theorem) Let ct = (a^^, 
of complex numbers such that for each i = 1, 
for each i 



s 



i 




(2.5) 



( 2 . 6 ) 



, Oj^) be any N-tuple 
N , I I = 1 , and 



(2.7) 



is a real number. Then, if (A^^, . . . , A^) is any stable point, it must 
be true that 



s^A^ + ... + Sj^ < 1 . ^ (2.8) 

(Note that since each s^ has modulus 1 and is real, each s^ must equal 
+1 or -1.) 

Proof. Let PC^ be the stable polynomial associated to (A^^, ..., 

Suppose a satisfies the hypotheses and define p(z) by 




p(z) = P(“iZ, . . . , a^z) . 



(2.9) 



Then p(z) is a one variable polynomial with real coefficients. Moreover 
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p(l) 



= 1 - 



K 

z 

i=l 



A^(a^, 




K 

A. s , . 
i 1 

i=l 




( 2 . 10 ) 



If 



K 



E 



SiAi 1 1 



( 2 . 11 ) 



then 



p(l) 1 0 . 



( 2 . 12 ) 



But p(0) = 1, and p(z), as a function of a real variable, is clearly 
continuous. Therefore by the intermediate value theorem for real functions, 
there must exist a real number t , 0 < t 1 , such that p(t) = 0. But 
then 

P(aj^t, ..., = 0 (2.13) 

and for each i = 1, ..., N , ^ 1 , a contradiction to the assumption 

that PCz) is stable. Therefore it must be true that 




.A. 

1 1 



< 



1 . 



(2.14) 



Theorem 2.3 (S3nnmetry Theorem) Let a = (a^, . . . , a^) be any N-tuple of 
complex numbers such that for each i = 1, ..., N , |a^| = 1 , and 



s . 

X 



s a 



(2.15) 
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is a real number. Then a point (A^, . . . , 

stable. (Note as before that each equals +1 or 

- 1 .) 

Proof. Let P(Zj^, ..., be the polynomial associated to (A^^, ..., • 

Define P* (z^^, ...» z^^) by 

P*(Zj^, ...) z^) = P(aj^z^, •.«> * (2.16) 

It is not difficult to check that P' is the polynomial associated to 

P'(6j^, .... = 0 iff P(8j^j ...» 8^) = 0 where 8j^ = and for 

each i |6!i = |ct^8^| = i3^|» since |a^| =1 . 



. .., Sj^Aj,) . Moreover P' is stable iff P is. For 



A^) is stable iff the point 



Theorem 2.4 (Sufficiency Theorem) Let (A^^, . . . , A^) be a point such 
that 




(2.17) 



Then (Aj^^ . . . ^ 

Proof. Consider the polynomial (1.7) 



A^) is stable. 



P(z^, 







If P(8) = 0 where |8^| ^1 > then 




(2.18) 



(2.19) 



Therefore 



1 = E A. < E l^il < 1 

i=l i=l 



( 2 . 20 ) 



a contradiction. 
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Corollary 2.1 Let A^,) be a point such that A^ ^ 0 for all 

X • Then ^A^^j • « • > 

K 

X) \ 1 • (2.21) 

i=l 

Proof. Apply Theorem 2.2 and Theorem 2.4. 

Theorem 2.4 and its corollary can be given simple geometric inter- 
pretations. The region of points (A^, ..., 

1Ai1 + ... + |a^1 < 1 (2.22) 

describe the K-dimensional 'diamond', centered at the origin and whose 
points lie along the axes. Theorem 2.4 states that the K-dimensional 
diamond is wholly enclosed by the region of stability. The Corollary 
states that in the positive 'quadrant', the region of stability always 
coincides with the 'diamond' (see Figure 1). In the following, the 
Symmetry Theorem will be applied to show that regions of stability also 
satisfy certain geometric symmetries. 

3. Symmetry 

As the proof of Theorem 2.3 shows, symmetries, (in this case sequences 
of reflections) between the points of the region of stability arise from 
transformations of the variables of the associated polynomials. The 
transformation of variable is given by 

(z^, . . . , • * • > ~ . . . , cijjZjj) (3. 1) 

where (a^, ..., a^) is a complex vector satisfying the conditions that 
for each i , |a^| = 1 and 




A^) is stable iff 



9 



a 



(3.2) 



s . = 



is a real nuinber. Such a stability invariant transformation transforms the 
coefficients (A^, 

always equals +1 or -1, this transformation is geometrically interpreted as 

a sequence of reflections of the axes. 

N 

At least 2 vectors a can be obtained by choosing each equal 

to +1 or -1. The symmetries of the coefficient space thus obtained, 
however, are not necessarily distinct, and may not include symmetries which 
can be obtained by allowing the to be complex. However, symmetries 

so obtained are more easily studied and for this reason Xi?ill be called 
simp le symmetries. In the following we will restrict our attention to the 
study of simple symmetries. 

The transformation (symmetry) of the coefficient space determined by 

the vector "a can be described completely by the vector s = (s^, ..., s^) 

X K 

specified by eqn. (3.2). Each coordinate of the vectors *a = (a^, ..., ot^) 
and ^ = (s- , ..., s^) equals +1 or -1. However for reasons which will 
become clear, it will be more convenient to use 1 in place of -1, and 0 
in place of 1. The multiplicative relation (eqn. 3.2) between the vectors 
's and ^ now becomes the following additive one in modulo 2 arithmetic. 






into 









Since each s . 



= "a • (mod 2) (3.3) 

where the * • ' represents the vector dot product of the 0, 1 vector "a 
with the integer vector index . If the N x K integer matrix J is 
defined by 

J = (t/. .... ij) (3.4) 
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then the vector s which describes the simple symmetry arising from the 
change of variables described by the vector a satisfies 

^ J (mod 2). (3.5) 

N 

Moreover the 2 possible choices for a can now be viewed as the 
elements of the N-dimensional vector space over the Galois field, GF(2), 
of two elements, a tool familiar in algebraic coding theory (Berlekamp 
[1968]). The simple symmetries can now be easily classified using the 
linear algebra of these vector spaces. 



Theorem 3.1 The set of K-dimensional vectors over GF(2) which describe 
the set of simple symmetries of the region of stability is the set of 
vectors spanned, modulo 2, by the row vectors of the matrix J . 

Proof. Each vector "s arises by eqn. 3.5 from an "a . Each 
a = (a^, ..., Oj^) can be written as 

a = a^e“j_ + ... + a^^-e^ (3.6) 

where 

SjL ~ C0> • • • > 1> • • • > 0) (3* 7) 

with a 1 in the i^^ coordinate. But then 



■^ = a^(^j^J) + ... + a^^Ce^^J) 



by linearity; and for each i , "e^J is the i^^ row of J . 



(Note also that by the above theorem the set of simple symmetries forms 
an abelian group.) 
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Corollary 3.1 The number of simple symmetries equals 2^ , where L is 
the modulo 2 rank of the matrix J . 

Proof. Immediate. 

Example 3.1 

Consider the filter (eqn. 1.8) studied by Huang. The vector indices 
in this case are 

= (10) 

^2 = (01) (3.8) 

= ( 11 ) . 

The matrix J is therefore 



J 



1 0 
0 1 




(3.9) 



2 

The modulo 2 rank of J is clearly 2. Thus there are 4 (=2 ) simple 
symmetries. Each symmetry is described by an element of the row space of 
J . They are: 



= (000) 

S2 = (101) 

i*3 = (011) 
i*4 = (110) 

which describes the symmetries 

(A^jA2>A2) (A^jA2>A2) 

(~Aj^ J A2 > “A^) 

(A2^)*A2j A^) 
(~Aj^>~A2> A^) * 



(3.10) 



(3.11) 
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Note that in Figure 1 , there are only two different basic shapes for the 
region of stability in the eight different quadrants of the coefficient 
space, and that consistent with the above symmetries each shape occurs 
symmetrically in four quadrants. 

Example 3.2 

Consider a filter whose transfer function is 






(3.12) 



where 



P(Zl»Z2» ^3*^4^ 





In this case the vector indices are 



= (1301) 
( 1221 ) 





(1331) 



(3.14) 



3^4 = (0331) 
T 5 = (1031) 



And the matrix J is therefore 



J 




0 2 3 3 : 



1110 ] 



3 2 3 3 C 



1111 ] 




(3.15) 



Calculating modulo 2, and row reducing 
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J 



(3.16) 




1 

0 

0 

1 



1 

1 

1 

1 



0 

1 

1 

1 




-> 




1 1 
1 0 
0 1 
0 0 



0 

1 

1 

1 




Therefore J has rank 4, and there are 2 = 16 possible distinct simple 

symmetries. Moreover each possible symmetry can be described by a modulo 
2 sum of the rows of J . For example adding every row we obtain 



■s = (1 0 0 1 1) 



(3.17) 



which corresponds to the s 3 nnmetry 



(Aj^» A2> Ag) (“A^> A 2 » A^j ~A^ , “A^) (3.18) 

of the coefficient space. Thus, for example, since we know that in the 
positive quadrant the shape of the region of stability is the part of the 
diamond in that quadrant it follows that in the quadrant which corresponds 
to the above symmetry, the region of stability is again the part of the 
diamond in that quadrant. 

Example 3.3 

Consider the class of polynomials above without the last term. 

= 1 - - k^z^z^z^z^ - A2Z^Z2Z2Z^ 

- A^z^z^z^ . (3.19) 

Since the rank of J will still be four, and the dimension of the 
coefficient space is four, it follows that every quadrant is symmetrical 
to every other by an appropriate simple symmetry. That is, every possible 
change of sign will occur among the simple symmetries. It follows that 
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the region of stability is the diamond, and 




< 1 



is a necessary and sufficient condition for the stability of polynomials 
of this type. 
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